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A LIOUVILLE-TYPE THEOREM FOR THE 3-DIMENSIONAL 
PARABOLIC GROSS-PITAEVSKII AND RELATED SYSTEMS 


QUOC HUNG PHAN AND PHILIPPE SOUPLET 

Abstract. We prove a Liouville-type theorem for semilinear parabolic systems of the 
form 

m 

dtUi AtAj — ^ Uj , % — 1,2,..., m 

j=i 

in the whole space x R. Very recently, Quittner [Math. Ann., DOI 10.1007/s00208- 
015-1219-7 (2015)] has established an optimal result for m = 2 in dimension N < 2, 
and partial results in higher dimensions in the range p < N /(TV — 2). By nontrivial 
modifications of the techniques of Gidas and Spruck and of Bidaut-Veron, we partially 
improve the results of Quittner in dimensions N > 3. In particular, our results solve 
the important case of the parabolic Gross-Pitaevskii system - i.e. the cubic case r = 1 
- in space dimension N = 3, for any symmetric (to, to)- matrix (fiij) with nonnegative 
entries, positive on the diagonal. By moving plane and monotonicity arguments, that 
we actually develop for more general cooperative systems, we then deduce a Liouville- 
type theorem in the half-space R;^ x R. As applications, we give results on universal 
singularity estimates, universal bounds for global solutions, and blow-up rate estimates 
for the corresponding initial value problem. 


1. Introduction 

we study the semilinear parabolic system of the form 

m 

— Au t — f3 i jU r l v r * 1 , i£fl, tel, i = l,2,...,m, (1) 

3 = 1 

where r > 0, Q is a domain of R iV , / is an interval of M, m > 2 is an integer and B = (/3^) 
is a real m x m symmetric matrix. We assume throughout, unless otherwise specified, 
that 


In this article, 

duj 

dt 


fiij > 0 for all i ^ j i and /3a > 0 for all i. (2) 

The system (jT]) can be used to describe heat propagation in a m-component combustible 
mixture [4], in this case Ui represent the temperatures of the interacting components. 
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In the special case r — 1, system ([I]) can be seen as a parabolic counterpart of the 
m-coupled nonlinear Schrodinger system 

m 

- Aui = ^2/3 ij u i \u j \ 2 , i — 1, 2, m. (3) 

3 =1 

The cubic system (|3l), also known as the Gross-Pitaevskii system, arises in mathematical 
models for various phenomena in physics, such as nonlinear optics and the Hartree-Fock 
theory for Bose-Einstein condensation (see e.g. [3j [9j [34], Tj). In nonlinear optics, the 
solution Ui stands for the f-tli component of the beam in Kerr-like photorefractive media 
(see e.g. [I])- The positive constant fin is for self-focusing in the i-th component of the 
beam. The coupling constant fiij ( i j) is the interaction between the i-th and the j-th 
components of the beam. In the theory of Bose-Einstein condensation (see [15]), Ui are 
the corresponding condensate amplitudes, fin and fi^ are the intraspecies and interspecies 
scattering lengths. The case fi^ > 0 means that the interactions of states |z) and | j) are 
attractive. 

System ([T]) has been recently studied in various mathematical directions such as: the 
local and global existence 121120], Holder regularity [7], symmetry property 1221 EH, blow¬ 
up behavior [18], and Liouville-type theorems HUGZHIH]. Our main goal in this paper is 
to prove Liouville-type theorems for the problem (JT]) and then to deduce their important 
applications on qualitative properties of solutions. 

We recall that Liouville-type theorems are statements about the nonexistence of solution 
in the entire space or in half-space. In recent years, the Liouvillc property has been refined 
considerably and has emerged as one of the most powerful tools in the study of initial 
and boundary value problems for nonlinear PDEs. It turns out that one can obtain from 
Liouville-type theorems a variety of results on qualitative properties of solutions such 
as: universal, pointwise, a priori estimates of local solutions; universal and singularity 
estimates; decay estimates; universal bound of global solutions, initial and final blow-up 
rates, etc..., see [2H 55] and references therein. In addition, it was shown in [35] that the 
parabolic system ([1]) can be used in the study of solutions of the corresponding elliptic 
problems, provided one can show suitable a priori bound of the global solutions of ([1]). 
This a priori bound property is a consequence of the Liouville-type theorems. 

Let us recall the elliptic counterpart 

m 

—A Ui = fijjU r i u r i +1 i seK 1 , i = 

3 = 1 

This system has attracted much attention of mathematicians in recent years, especially 
for the cubic case r = 1, see e.g. (3oui3i m\ edi is nzj for more references. Concerning 
the Liouvillc property, it is well known that the Liouville-type result for (|4]) plays an 
important role in the study elliptic problems as well. The optimal Liouville-type theorem 
for nonnegative solutions of (QJ was completely proved by Reichcl and Zou [32] (see also 
[14]) via moving sphere techniques, under the optimal Sobolev subcritical range p < ps, 


( 4 ) 


1 dui 

yf-l dt 
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where 


and 


p := 2 r + 1 


Vs '■= 


N+2 

N-2 

OO 


if N >3, 
if AT = 1,2. 


We note that, if we remove the positivity assumption on the diagonal of B , then problem 
(Hj) . and hence ©, may have many semi-trivial solutions due to system collapsing, i.e., 
solutions with one or more components being zero. For example, if &\ \ = 0 then U = 
(C, 0,..., 0) is a solution of problem (0}. And indeed, the study of the system (0]) becomes 
more delicate due to the existence of semi-trivial solutions, especially in the applications 
of Liouville-type theorems (cf. pa, USEE] and see Remark 00] below). 

For the corresponding parabolic problem (0]), the Liouvillc property is much less under¬ 
stood. Recall that, even in the scalar case, i.e. for the classical nonlinear heat equation 


u t — A u = u p , 


(5) 


it is by now not completely settled. Indeed, the Liouville property for (0]) is conjectured 
to be true under the optimal condition 1 < p < p$ (and this is known to hold in the class 
of radiallly symmetric solutions [251 [23] ). but this has been proved so far only under the 
stronger restriction 1 < p < ps where 


Pb(N) : = 


N(N+2) 

(JV-1)2 

OO 


if N > 2, 
if N = 1, 


( 6 ) 


or, very recently EZ|. for N = 2. One of the main difficulties is that the techniques of 
moving planes or moving spheres do not work as in the elliptic case [32]. As for system (Jl]) 
under assumption (J2J), only some partial cases are known: 

• First, one can easily obtain a scalar parabolic inequality d t z — A z > Cz p for 

z = an d deduce the Fujita-type result of problem (0J, namely there is no 

nontrivial nonnegative solution in x K + ifl<p<l + ^. 

• In the case m = 2, the Liouville-type theorem for (0]) has been recently proved 
in [2T| in dimension N = 1 and for radial solutions in any dimension if p < p$- 
More recently, Quittner [27] has proved the optimal Liouville-type theorem in 
dimensions N < 2, and has also given a partial result in dimension N > 3 under 
the condition p < jfh- r The main tools in [27] are scaling argument and energy 
estimates. 

• Under the assumption m = 2, f3a = 0 for i = 1,2 and /3 i 2 > 0, the Liouville- 
type theorem for positive solutions of problem (0]) can be shown via comparison 
technique (see [28] and [31] 09] for the elliptic case). More precisely, by taking 
the difference of the two equations and suitably using the maximum principle, we 
may show that u — v and thus reduce the system to a scalar equation. 






4 


Q. H. PHAN AND PH. SOUPLET 


In this paper, we shall use a different approach to establish a Liouvillc-type theorem 
for problem (OQ) in the whole space in a larger range of p and for any m. We shall then 
treat Liouville-type theorems in the half-space by reduction to the whole space case. 

2. Liouville type results 

Our main result in the whole space case is the following. 

Theorem 2.1. Let m > 2, r > 0 and assume p := 2r + 1 < ps(N), where ps(N) 
is defined in ®- Let B satisfy (0|). Then system m has no nontrivial, nonnegative 
classical solution in x R. 

Remark 2.1. a) We stress that p B (fV) > when N > 3, and our result is a partial 
improvement of Quittner [27j in higher dimensions. In particular, it solves the important 
case of the parabolic Gross-Pitaevskii (cubic) system where r = 1 and N = 3. 

b) If N = 2 then, by [27] . the conclusion of Theorem \2.1\ is actually true for all r > 0 
(the result in [27] is formulated only in the case m = 2 but the proof is valid for any 
m > 2). 

c) Our proof of Theorem \2. 1\ relies on nontrivial modifications of the technique developed 
by Bidaut- Veron [5] for the scalar nonlinear heat equation. The latter was an adaptation 
of the celebrated method of Gidas and Spruck [12] for elliptic equations (see also 0 for 
some particular elliptic systems). It is based on nonlinear integral estimates and Bochner 
formula. This technique is completely different from that of [27] , which relies on scaling 
and energy arguments. 

d) In [TS], Merle and Zaag proved Liouville-type theorems for the so-called ancient 
solutions of the system 

U t — AU = F(\U\)U, (7) 

with F(\U\) ~ |I7| P_1 as \U\ —* oo, under the assumption p < ps and 

Ui(x,t) <C(T-t)~ l/(p - l) . ( 8 ) 

Namely they showed that any solution of £/]j in ~R N x (—oo,T) which satisfies is inde¬ 
pendent of the space variable. In the special case F(\U\) = \U\ 2 (p = 3), Proposition \3.1\ 
below (which is a consequence of Theorem \2.1\) guarantees that §3) holds if 3 < ps(N). 
Thus, the estimate m is always true if N < 3. 

We now turn to the case of a half-space R^ = {i£ M. N : Xi >0}. 

Theorem 2.2. Let r > 1 and assume either N < 3 or N = 4 and p = 2r + 1 < Pb(3) = 
15/4. Let B satisfy (d[). Then the problem 

{ du m 

- Aui = y ^ (3 ij u r i u r j +1 , x G R+, te R, i = l,2,...,m, 

J=1 

Ui =0, x G OR”, tE R, i = l,...,m, 

has no nontrivial, nonnegative, bounded classical solution. 


( 9 ) 
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Remark 2.2. The upper restrictions on N in Theorem \2.2\ are consequences of the con¬ 
dition r > 1, which is required by the use of the moving plane method in the proof of 
Theorem 1 2. ‘A 


Theorem 12.21 is a consequence of the following result which can be applied for more 
general cooperative parabolic systems. 


Theorem 2.3. Consider the following system 

Q Ui 

-jC-- Aui = fi(ui, ...,it m ), X e R+, t e R, i = 

Ui = 0, x G teK, i — 1, ..., m. 

Assume that fi : [0, oo) m —$■ M are C 1 -functions satisfying: 

(Ri) /i(0,...,0) =0, i — 1, ..., m, 

dfi 

(H 2 ) -—A(ui,...,u m )>0, for all [u \,..., u m ) G [0, oo) m andalliy^j, 


( 10 ) 


du 


dfi 


(H 3 ) - °’ f° r al1 

j =i ° j 


and that 

(Hf) any nontrivial, nonnegative, bounded solution of fTW is positive in x M. 
Then any nontrivial, nonnegative, bounded solution U = (uf) of fflft) is increasing in X\: 


duj 

dxi 


(x, t) > 0, x e tel, i — 1,..., m. 


( 11 ) 


Theorem 12.31 is an analogue for systems of ( 25l Theorem 2.4 (cl)] for scalar equations. 
The proof follows the idea in [25] which is based on a moving plane technique. However, 
significant additional difficulties arise in the case of systems. This leads to the introduction 
of the assumption (Hf) (which, in turn, is necessary for the conclusion (flTj) to hold). Then, 
in order to deduce Theorem 12.21 from Theorem 12.31 we use an induction argument on the 
number m of components. 

The outline of the rest of the paper is as follows. In section 3, we give applications of 
our Liouville-type theorems, namely universal singularity estimates, including initial and 
final blowup estimates, as well as universal bounds for global solutions. Section 4 is then 
devoted to the proof of the Liouville-type Theorem 12.11 in the whole space, and section 5 
to the proofs of Theorems 12.21 and 12.31 in a half-space. Finally, a version of the maximum 
principle for cooperative systems, suitable to our needs, is given in Appendix. 


3. Applications of Liouville-type results 


As a first application of Theorem 12.11 we obtain universal singularity estimates in time 
and space, including universal initial and final blowup estimates in the case H = W N . 
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Proposition 3.1. Let r > 0 and assume p := 2r + l < ps(N) if N > 3. Let B satisfy (0|). 
There exists a universal constant C = C(N,p, B) > 0 such that, for any domain 12 ofR N 
and any nonnegative classical solution U of fj\) in Ll x (0 ,T), there holds 

Ui {x , t)<C (V 1 ^- 1 ) + (T — t)- 1/(p ~ 1} + dist“ 2/(p_1) (x, <912)) , (12) 

for all (x, t) G 12 x (0, T) and i — 1,..., m. In particular, in the case 12 = , we have 

Ui (x, t)<C (r 1/(p_1) + (T - t)-V(p-i)) . (13) 

Note that Proposition 13.II covers in particular the Gross-Pitaevskii case r = 1 in dimen¬ 
sion N = 3. The proof is based on a reduction to Theorem 12.II bv rescaling and doubling 
arguments. Since it is completely similar to that in [25l Theorem 3.1] for the scalar case, 
it is therefore omitted. We stress that the proof only requires Theorem 12.11 for bounded 
solutions (and Theorem 12.11 in the general case is finally obtained as a consequence of 
Proposition 13. T)) . 

Remark 3.1. (a) It is clear that if the positivity condition of the diagonal of B is removed, 
then the universal singularity estimate TT2 1) fails due to the existence of arbitrary large 
semi-trivial solutions. Indeed, if [In = 0, then for any A > 0, the constant function 
U = {A, 0, ...,0) is a semi-trivial solution of problem (ED in any domain 12 x (0, T). 

(b) In the elliptic case, it is sometimes possible to prove universal estimates of positive 
solutions in spite of the existence of arbitrary large semitrivial solutions. This is for 
instance the case for Dirichlet problems associated with system |^) when r — 1, m — 2, 
[In = /?22 = 0, /?i 2 > 0 and N < 3 (see [H5J Theorems 1.1 and 6.1]J. However the following 
counter-example shows that even this fails in the parabolic case: 

Let r > 1 and assume that fin = 0. Let A > 0 be fixed. For any e G (0, A), we denote U £ 
the maximal classical solution of system m in B i x [0, T) with initial and boundary value 
U e = (A,e, ...,£). Then u ei > e > 0 by the maximum principle. Also, U 0 = (A, 0, ...,0) is 
a global solution of system (EP in B x x [0, oo) with initial and boundary value (A, 0,..., 0). 
Since the nonlinearity is Lipschitz (r > 1), it follows from the continuous dependence of 
solutions with respect to initial and boundary data that, for e G (0, e 0 (A)) with £q(A) > 0 
small enough, the solution U £ exists in B 1 x [0,1]. Since A > 0 is arbitrarily large, the 
universal estimate / EH) thus fails also for positive solutions. 

For the next application of our Liouville-type theorems, we consider the initial-boundary 
value problem: 

-^-Aui + XiUi = y x G 12, f G (0,T), i = l,...,m, 

< j=1 (14) 

Ui =0, x G <912, f G (0,T), i = l,...,m, 

Ui(x, 0) = Mo,*, x G 12, i — 

where Uq = («o,i) are non-negative functions in (Go(12)) m , Aj G 1 are constants, and 12 
is a regular domain (possibly unbounded) of M A . It is well known (see, e.g., Amann [2]) 
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that problem (fT4l) has a unique (mild) solution 

tf(s)ec([o,r);(c„(n)r) 

with maximal existence time T = T max (f/o), and U is a classical solution for t G (0, T). We 
have the following universal estimates for global solutions, as well as universal initial and 
final time blow-up rates. Again, this covers in particular the case r = 1 in dimension N = 
3. 

Proposition 3.2. Let r > 1 and assume either N < 2 or N = 3 and p = 2r + 1 < 
Pb( 3) = 15/4. Let B satisfy (U|). Let U be a nonnegative solution of [Tf\ ) in D x (0, T). 

(i) If T < oo then there holds 

Ui(x, t) < C (l + _|_ (j 1 _ ^-VIp- 1 ) j ^ x G 0 < t < T, i = 1, ..., m, 

where C = C(Q,p,B). 

(ii) If U is global then there holds 

ufx , £) < C (l + 5 x G t > 0, i = 1, ..., m, 

where C = C(Q,p, B). 

The proof of Proposition 13.21 is completely similar to that of |25] Theorem 4.1], based 
on rescaling and doubling arguments, and Liouvillc-type theorems. Namely, we use the 
Liouvillc-type Theorems 12.II in the whole space, and Theorem 12.21 in a half-space. 

Remark 3.2. The upper restriction on N in Proposition \3. 6 L\ is a consequence of the 
condition r >1, required by Theorem \2. 6 A However, if we consider only radial solutions 
in a symmetric domain 11 (i.e. the whole space WL N , a ball, an annulus, or the complement 
of a ball), then Proposition \3.2I is true under the weaker assumptions r > 0, N > 1 and 
p < ps (see [21, Section 2}). 


4. Proof of Theorem 12.11 

For the sake of simplicity, we denote by J the integral f B dxdt. The key step is the 
following Lemma. 

Lemma 4.1. Assume that p < ps(N) and B satisfies (|I|). Let 0 < p G T>(Bi x (—1,1)) 
and U be a positive classical solution of m in B\ x (—1,1). Denote 

P m 

Ii= I puf |Vu,| 4 , / = £>, 

1=1 

nn 2 m 

^(j]ftxy +1 ), i = E £ <- 

j =1 i =1 
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Then there holds 


m r> m r, 

I + L<cJ2 / (p{\dtu i \uZ 1 \Vu i \ 2 + \d t u i \ 2> j +CJ2 / |Avj||Vui| 

i= 1 ^ i= 1 ^ 

m p m 

+ Cj^ / (E^r+IW+^lV^lVfp.Vuil 

i=l ^ .7—1 


+c 5 ^ / ixftx 

*j=i J 

where C = C(N,p, B ). 

Proof. Step 1. Preparations. Denote 


r+l-.r+l 


(15) 


/ m 

Ai/j, J ^ ^ 

i=i 

/ m 

p{ An,) 2 , K = Y,Ki. 

i =1 


Applying [29. Lemma 8.9] with q = 0, — 1 7 ^ k < 0, we have 

A; + ijfcA + 


N -l 
N 


N 


N 


< - j \X7ui\ 2 Aip + J (A-Uj — ku t 1 |V?ii| 2 ) Vwj.V^. 

(We stress that this is true if Ui is any positive C 2,1 function, with no reference to the 
PDE system (JTj) at this point.) Therefore, 


N - 1, \ N + 2, IV- 1 

k + l)kl H-77 —kJ -77—A 


N 


N 


N 


m , 1 

/ f-|Vnj| 2 A(/? + (An* — Vuj| 2 )Vni.V</9 

i=l ^ '■ 


(16) 


Step 2. Estimate of J and K. We claim that 


1 m « m 

-•x-£ M£ax 

^ *=i J j=i 


, 2 

r ^ +1 ) + 


2 p(r + 1 ) 


<Pt 


v ftx 


r+l,.r+l 
U 3 


i,j =1 


„ m „ m 

^ ^ ij = l ^ i=l 


m n. m 2 m /» 1 m /* 

k =J2j + J 2Jv\ d t u i\ 2 + —^[Yl I 

1=1 j= 1 2=1 hj = 1 J 


x + 1 y +I - 


(17) 


(18) 


and 
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Let us first establish (ITT]) . Using integration by parts, we have 


<p\Vu i \ 2 u r i - 1 u r j +1 = I pu r j +1 'Vu i .'V ( — 


= — - / ipulu^ 1 Alii - - I u r i u r j ¥1 \/ip.'S/u i 


i “j 


r + 1 


< pu r i u r jVu i .Vuj. 


(19) 


Using the Young inequality 2 u r i u , j 'Vu i .'Vuj < |Vuj| 2 w[ 1 u r j + | Vu^u^u 1 ^ 1 , it follows 
from (1T9|) that 

J <p\Vu i \ 2 u r i - 1 u r j +1 >~\J <pu r iUr +1 Aui - i j 


r + 1 
2 r 


^(iVUilX x ^ +1 + |VM i | 2 < +1 uJ x ). 


Consequently, 


rri ,, m ^ ra „ 

E Pa / ^iv« i | 2 <-x +i >~J2 --J2 

2j=l ** i,j= 1 ^ 2,j=l ^ 

m p 

■ hi / ^(|Vui|\ r_1 ^ +1 + |Vu/i^+yp 1 ), 

j,j=l J 


r + 1 


hence 


m r, ^ m p ^ m p 

£ Pij / ¥>|Vw i | 2 <- 1 ^ +1 >--J2 ^9<«r lAu >--E / /3y<«- + 1 Vv>-V«i 

2,7 = 1 ^ r 2,7=1 ^ r 2,j = l ^ 

m p 

^Pa / 

i,i=i J 


r + 1 




owing to the symmetry property /3„- = j3ji. Consequently, recalling p = 2r + 1, we obtain 


m p p m ^ p m 

Y1 Pij / ^|VM i | 2 <“ 1 u5 +1 >-- J2 <PPijUiUj +1 AUi -- J2 

2,,7 = 1 ^ y J i,j—l ^ ^ i,j=\ 


A J <«; + 1 v^.vu i . 


( 20 ) 


Now substituting the PDE (JT]) , written as 

m 

- Aui = ^2 Pij u i u ^ +1 ~ dtUi, 
j =i 


(21) 






10 


Q. H. PHAN AND PH. SOUPLET 


in the definition of J, it follows from (T20l) that 

m r> 171 

-j =y / (p\'vui\ 2 /3i j u r i ~ i u r j +i - 

M =1 ' i =1 ^ 

“ 1 m r 

;-a«o --'521 ^ u >? 1 iv^.v^i 

P i,7=l 


r M r+! 


-| I/O n HO 

£ -£ / KEftx 

^ t=l J J=1 

m « 

~52 v\dtUi\u~ l \\7u 
i=1 ^ 


and then, that 


r r+1 

i-«iV 


2 -- 

pJ 


■ l3 UiU^ +1 d t Ui 


, HO n HO 

F i=1 J j=1 

I 77L /» i/6 a 

-~52 Piju>r +i \v<p.vu i \ -J2 / ^i^iwr'ivMii 2 . 

p »,j=i ■' i=i J 

Next observe that, due to = /3ji , we have 

m m 

(r + 1) ^ fi i3 u\u r + x d t Ui + (r + 1) ^ ^u^u^dtu 


a.(E« 

»d=i 


*X + X +1 ' " 


*J=i 

m 


i,j=1 


2(r + !) 

»j=i 


Combining (|22|) and (1551) and integrating by parts in t, we obtain 


m « m 

' J -- v 52 / 

^ i=i j j= i 


\ 2 i 

« +i i + 


' 7 2p(r + 1) 


Vi 


V ftx 


r+1 r+1 
U 3 


»J =1 


m r. m r. 

52 / - E / 

i,j=l ^ i=1 ^ 


i.e., estimate mm. 

Next consider K. By substituting (l2Tj) in the definition of fT, we have 

k = 51 / + 52 <^| 2 - 2 £ / yPi j u r i u r + 1 d t u i . 

i=1 J j=1 i=1 ^ ij=l ^ 

Using (|23|) again and integrating by parts in t, we obtain (1151) . 
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Step 3. Conclusion. It follows from (1171) and (1181) that, for any — 1 k < 0, 

2 


iV + 2, T TV-1 

-kJ -——A > 


N 


N 




r r+1 


lib n lib 

Y \ t{Y^a 

i= 1 J j =1 

m r, m « 

-CE ) wftx + S r+1 - c E / 

2,7 = 1 °' 2,7 = 1 "' 


>»J= 

m 


m 


C Y ^\d t Ui\u~ l \S7ui\ 2 - / ^|< 9 t «/ 2 . 

i=i *=i •' 


(24) 


Since p < pb(N), we can take k > —N/(N — 1) close to — iV/(iV — 1) such that 


(V* 11 ) 1- ^ 1 and y^r ( “ fc) - > °- 

(If iV = 1 then we take any —1 7 ^ k < 0.) The lemma follows from (1T61) and 


N — 1 


□ 


Lemma 4.2. Assume that p < ps(N) and B satisfies (0J). Let U be a positive classical 
solution of m in B 1 x (—1,1). Then 



,1/2 


'-V 2 “T 


A+ 1 ") 


1=1 


dxdt < C(N,p, B ). 


(25) 


Proof. We follow the argument as in the proof of [2H1 Proposition 21.5]. One can choose 
the test-function </? such that <p = 1 in Ph /2 x (—1/2,1/2), 0 < ip < 1 and 

|V^| < CV (3p+1)/4p , |A</?| < Cp {p+1)/2p , \p t \ < Cp {3p+1)/4p < Cp {p+1)/2p . (26) 

Recall the notation J = f B f\ 1 dxdt. By the proof of [29| Proposition 21.5] (see formulae 
(21.10) and (21.11)), for any s > 0 and any positive function u G C 2,1 (i?i x (—1,1)) , we 
have 

J\Vu\ 2 (\Ap\ + + \p t \) < £ J p(u~~\Vu\ 4 + u 2p ) + C{e) 

and 

J [p\u t \u- l \\7u\ 2 + ( u p + |it*| + m _1 |Vu| 2 ) |Vu.V<^| + u p+1 \p t \) 

<£ J p(yC 2 \X7u\ l + u 2p ) + C(e) ^1 + J p\u t \ 2 ^j . 

Set /3 = min fin > 0. Using 

l<2<m 


and 


m „ m r, 

I + L>J2 fim 2 |Vu,| 4 + /? 2 E / 

i =1 ^ i =1 •' 




2 p 


«>r <«?+«?, 
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it follows that 


V?/i|“(|A<^| + ip 1 1Vv?|^ + \^pt\) A s(I + L) + C(e), 

m 

(p\d t u i \u~'-\Vu i \ 2 + (E /%« +1 + \dtUi\ +M“'|VM i | 2 )|Vu l .V^| 

3 = 1 

< £(/+ L) + C(e) ^1 + J ip\dtUi\ 2 ^, 

m „ 

(E^< + 1 -r> t |) < E + £) + C^e )(i + / ¥>|&«i| 2 )- 


(27) 


i=i 


On the other hand, using (1231) and integrating by parts in space and in time, we have 


m r, m p m 

E / v\dtUi \ 2 = E / + E Ad 

i=l J \ * j j ij =1 

i r 


= / '•Pt 


E 

%— 1 


|V?/i 


< C 


2(r + 1) 

/ /» l Tfl 

n(ei v “‘I 2+ e« 

i=l i, 7=1 


<Pt 


E ax 

*d=i 


r+1 r+1 


r+1 r+1 
a 3 


110 p 

- E / (9tUi)V<p.Vui 

i= 1 ^ 
m ,, 

E / (dp+V^.Vu* 

i=i ■' 


r+i r+i 

lj LL J 


m p .. m r. 

E / 7X M i| 2 + a E/ / < 7’“ 1 l V ^| 2 |V*h| 

i=l •* i= 1 ^ 


(Note that, since n* is positive, it is smooth enough so that the above calculations are 
justified.) Therefore, 


m p m p m 

E / xx 2 < c 'E / n(iv+ 2 +E/% 

*=1 ^ i=l ^ j=l 




E f + 1 |VH 2 |V«, 

i=l ^ 



+ C(iV,p, B)e(I + T) + C(s). 


(28) 
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By (1261) and Young’s inequality, 

j M< + \ r+1 < 2 eJip(u iUj ) p + C(e) J<p-( p+1 ^ p -V\<p t \ 2p ^ p -V 

<e J p(u 2p + u 2p ^ + C(e) |„-(^ 1 )/(^i)|^/(p-x) 

- + oF^' + ^( £ )- 

r'ii i jj 

Hence, 

J MAX +1,u i +1 — A? (^2 + ^r) (-^ + -k) + C( 5 )- ( 2 9) 

Combing (12g|) and (J22J), we obtain 

m p 

/ ‘f\d t u i \ 2 <C(N,p,B)s(I + L) + C(e). (30) 

i=i ^ 

Therefore, it follows from (T15T) . (127|) and (1301) that 

I + L < C{e) + C(N,p, B)e(I + L ). 

By choosing e sufficiently small, we obtain I, L < C and the Lemma follows. □ 


Proof of Theorem \2.1[ We first consider the case of bounded solutions. Assume for 
contradiction that U = [uf) is a nontrivial, bounded, nonnegative solution of (P in 
~R n x M. For each i, since the component u* is a supersolution of the heat equation, it 
follows from the strong maximum principle that either iq is positive in W N x R, or there 
exists t 0 G M such that — 0 in M, N x (—oo,£ 0 ]. In the latter case, since U is bounded, 
we have d t Ui — Aw, < Cui for some constant C > 0 and the maximum principle then 
guarantees that tq = 0 in x (t Q , oo), hence iq = 0. By removing all the components 
which are identically zero and relabeling, we may assume without loss of generality that 
Ui > 0 for i = 1 ,..., m (with m > 2 , since nonexistence in the scalar case is already known 
by [5], see also |29] Theorem 26.8]). 

Now, for any R > 0, we rescale 

Vi(x,t) = R 2 ^~ r> u i {Rx,R 2 t). 


Then V = (vf) is also a solution of (JT|). Since all Vi are positive, we may apply Lemma 14.21 
to deduce 


m 


'\y\<R/2 


' \s\<R 2 /2 “ 


m 

E 

3 = 1 


axy 1 ) 


2 

(y, s)dyds 


r. r. mm 

/ / £(2> 

J\x\<l/2 J\t\<l/2 l=1 \ / = 1 


— jj>JV+2-4p/(p-l) 

< CR N+2 ~ Ap ^ p ~ 1 ^ 


XX 1 ) 


2 

(x, t)dxdt 
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Since p < ps < ps, by letting R —> oo we get Pij u i u1 j +1 ) 2 = 0- Since /3a and 

Ui are positive, this is a contradiction. This proves Theorem 12.11 in the case of bounded 
solutions. 

Finally, to treat the general case, we recall that the Liouville-type property of The¬ 
orem 12.11 for bounded solutions is sufficient for the proof of Proposition 13.11 (see the 
paragraph after Proposition 13.ip . But after a time shift, formula (fT31) in Proposition 13.11 
guarantees that any solution of ([I]) in R N x (—T, T) has to satisfy u(x, t) < CT -1 /(p-i) in 
x (—T/2, T/2). The conclusion then follows by letting T —* oo. □ 


5. Proof of Theorems 12.21 and 12.31 

Proof of Theorem 13.31 It is done in three steps. 

Step 1. Notation. For A > 0, let 

T A = {i 6 M a : 0 < X\ < A}. 

Denote 

V\Ui(x , t) := Ui( 2A - xi,x', t ) - Ui(xi,x', t ), 
where x 1 = (# 2 , ■ Let V{ = V\Ui, then V = (vi) satisfies 

m 

d t Vi - A Vi = ^2 c ij v v t) e T A x M, % = 1,2,..., 

3 = 1 

Vi = 0, X\ = A, x' e M^ -1 , t e R, i = 1,..., m, 

Vi >0, X\ — 0, x' G M^ -1 , t e R, i = 1,..., m, 


m. 


where 


c A - = 

*3 


" l dh_ 

, dui 


(■U + sV)ds. 


We shall show that, for any A > 0, 

V > 0 for all (x, t) e T\ x M. 

For any positive g, A satisfying A^/g < 7r, we define the function 


which satisfies 


= sin - (-7T + y/q(2xi - A) 

—A h = qh, 

x e t a , 

h(x) > rj > 0, 

x e t a , 

|V/i| < y/q, 

x e T a , 


where rj = sin[|(7r — Ay / g)] E 


(31) 


(32) 

(S) 

(33) 


(34) 


1 We note that this choice of h is simpler than that in mm, owing to the different form of maximum 

principle used in the subsequent steps (which does not require W —> 0 at space infinity in (13^1) or (1451) 1. 
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Step 2. Proof of (S) for small A. Due to the boundedness of U, the coefficients cA are 
bounded (uniformly in A). So we may fix 7, q > 0 such that 


q > 7 + max sup . E 

1 A>o j =1 


(35) 


For any A > 0 small such that A 2 i r 2 > q, let h be given by (133]) and denote Wi := 
e^Vi/h. Then W = (wf) satisfies 

' m 

d t Wi - A Wi - ^.Wwi = (7 + c\ - q)Wi + ^ c^Wj, (.x , t) G T A x M, i = 1, m, 

3 = 1 

uij > 0, x G <9T A , t G K, i = 1, m. 

(36) 

By assumption (i/2), we see that cA is nonnegative for all i 7- j, so that the system (ITT]) 
is cooperative. Moreover, by the definition of q in (1351) . we have 

(7 + 4 - g) + 53 4 < 0, for all i = 1,m. 

We may thus apply the maximum principle for cooperative parabolic systems (see Propo¬ 
sition |(i. ll ii) in Appendix). The latter, applied to —W on T\ x (f 0 ,£) for any t 0 < t, 
guarantees that 

max sup w~(x, t ) < max sup wf(x, t 0 ), (37) 


where 2; := — min(2, 0). Consequently, 

< e 


max sup V -±SffA < e - 7 h-to) max sup !lM! < e - 7 (t-t 0 )MU 
* xgt a h(x) i o;gt a h(x) 77 


(38) 


Letting f 0 —> —00, we obtain L > 0 in Ta x M. Therefore, (S) holds for any A > 0 
sufficiently small. 

Step 3. Proof of (S) for large A. Next, we denote 

Ao = sup{/i > 0 : (S) holds for all A G (0,/i)}. (39) 

The previous argument shows that Ao > 0, and we shall prove by contradiction that 
Ao = 00. Assume Ao < 00. Then there is a sequence A*, > Ao such that X k —>■ Ao and the 
set 

Z k = {(x,t) G T\ k x M : there exists i such that V\ k Ui(x,t ) < 0} 
is nonempty. Set 

m k := sup max < maxuj(x, t), maxuj(2Afe — xi, x', t) >. 

(x,t)ez k 1 J 

We have the following two possibilities: 

Case 1: rrik > £0 for some e 0 > 0, up to a subsequence; 
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Case 2: m k —> 0. 

First consider Case 1. Then there exist sequences x\ G (0, A*,), z k G M A_1 , t k G R such 
that 

(x k , z k , t k ) G Z k and max | ma xui(x k , z k , t k ), max«j(2A fc — x k , z k , t k ) | > e 0 for any k. 

By extracting a subsequence of k, we may assume that there exists io G {1, such 

that, for any k, 

V Xk u io (x k ,z k ,t k )< 0 (40) 

and 

max | maxui(x k , z k , t k ), maxui(2X k — x k , z k , t k ) j > e 0 . (41) 

We may also assume that x k —* a for some a G [0, Ao]. Let 

u k (x, t) = Ui(xi, x' + z k , t + t k ), x = (xi, x') G R+, t G R, z = l, 

Since the sequence U k = (u k ) is uniformly bounded, using standard parabolic estimates, 
it follows that U k converges (up to a subsequence) in Cf^!(R+ x R) to a nonnegative 
solution U = ( Ui ) of ([TO]) . 

The definition of Ao implies that V Xo u k > 0 in T Ao x R for any i = 1 Hence, 

V Xo Ui > 0 in T Ao x R for any i = 1 Let v t := V Xo ili. Then (hj) is a nonnegative 

solution of the system 

/ m 

d t Vi - Avi = 'Y^c^Vj, (x,t) e T Ao x R, f = 1, 2,...,m, 

< j=1 

Vi = 0, X\ — A 0 , x’ G R^ -1 , t G R, i — 1, ..., m, 

^Vi > 0, x\ — 0, x 1 G R^ -1 , t G R, i — 1, ..., m. 

From gnD and (HI]) , we deduce that 

I/ Ao fq o (a,0,0) < 0 (42) 


maxMj(a, 0, 0) > £o or maxfh(2Ao — a, 0, 0) > Eo. (43) 

i i 

Owing to assumption (i/4), (jUfj) guarantees that 

Ui 0 (x,t) > 0, for all (x,t) G R+ x R, 
hence V Xo u io (0,x',t) > 0 for any x',t. Since z := V Xo u io satisfies 

( d t z - A2 > cf° iQ z, (x, t ) G T Ao x R, 

J z = 0, X\ — A 0 , x' G R^ -1 , t G R, 

[z >0, x\ — 0, x' G R^ -1 , t G R, 


we deduce from the (scalar) maximum principle that V Xo Ui 0 > 0 in T Ao x R. It thus follows 
from (]42]i that a = Ao- By the Hopf boundary principle, 

^dx^UiQ (Aq, 0,0) dxi Lao^jq (*^1 j 0) 0) |xi=Ao &xiz ( Ao> 0,0) > 0. 
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Consequently, d Xl Ui 0 (x i,0,0) is bounded below by a positive constant on an interval 
around A 0 and this remains valid if u io is replaced by u\ Q ~ u io . That is, there is S > 0 
such that 

d Xl u io (x l7 z k ,t k ) = > 0, xi e [A 0 -5, A 0 + 5] (44) 

for k sufficiently large. However, since 2A& — x\ > x\ both belong to [Ao — S, Ao + 5] for 
large k, (144)) contradicts V\ k Ui 0 {x\, z k ,t k ) < 0. 

We next turn to Case 2. We consider the system (j3Tj) for large k. We are going to apply 
the maximum principle again, this time taking advantage of the fact that the parabolic 
inequalities need be satisfied only on the possible positivity set Zk (see Proposition 16.11 ii) 
in Appendix). Recalling V = ( V\ k Ui ), it follows from rrik —> 0 that 

sup (\U(x,t)\ + \V(x,t)\) ->0 as k — y oo. 

(■ x,t)ez k 


In view of assumption (Hf) and the definition ([32]) . we deduce that 

m 

lirrisup < 0, where qk : = max sup 
k->oo ( x,t)ez k - =1 


Fix a large k and contants <7,7 > 0 such that q = qk + 7 < A fc 2 n 2 (« A 0 2 7 t 2 ). Like in 
Step 2, we consider Wi := e^Vi/h with h given by (1331) for A = A k- Let 

Di \= {{x,t) e T Xk X R; Wi(x,t) < 0} 


and note that C Zk- Then [wi) satishes 

' m 

d t Wi - A Wi - ^.ViUi = (7 + - q)wi + ^ % k wj, (x, t ) e A, i = 1, -,m, 

Wi > 0, (x, t ) e x R, i — 1,..., m. 


3 = 1 
3 & 


(45) 


By the maximum principle in Proposition 16.ll nh applied to —W, we obtain again (]37j) - 
(]38|) . and conclude that V > 0 in T\ k x R. But this is a contradiction with the nonemp¬ 
tyness of Zk- 

We have thus reached a contradiction in both cases, which proves that A 0 = 00 i.e., 
(S) holds for any A > 0. The Hopf boundary principle then gives 

2<9 3 , 1 'Uj((C, f) | a;i_=A 9xiV\ 0 Ui(^X,t^\ xl =\ 0 

for any A > 0 and any i. The theorem is proved. □ 


Proof of Theorem \2.2 l It is done by induction on m. For m — 1, it is reduced to the 
Liouvillc-type theorem for the equation ut — A u = u p in the half-space R^ x R. This 
was proved in [251 Theorem 2.1] for N < 2, or N > 3 and p < ps(N — 1). The result of 
Quittner 12?], in conjunction with [25l Theorem 2.4], guarantees the case N = 3. 
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Now fix m > 2 and assume that the theorem holds for m — 1. We shall prove by 
contradiction that the theorem then holds for m. Thus suppose that U = (iq,... ,u m ) is 
a nontrivial, nonnegative, bounded classical solution of ([2]). 

First we claim that 

any component w* is either identically zero or positive in R^ x R. (46) 

Indeed since tq is a supersolution to the heat equation, it follows from the strong maximum 
principle that either tq is positive in R^ x R, or there exists to £ R such that tq = 0 in 
R+ x (—oo, t 0 ]. In the latter case, since U is bounded, we have d t Ui — Atq < Chq for some 
constant C > 0 and the maximum principle then guarantees that tq = 0 in R^ x (t 0 , oo), 
hence tq = 0. 

Next we claim that actually 

Ui > 0 for allz = 1, ...,m. (47) 

indeed, otherwise, by fl46l) . there is a component tq 0 which is identically zero. By removing 
this component, we then obtain a nontrivial, nonnegative, bounded classical solution of 
problem ([9]) for (m — 1) components. But this contradicts the induction assumption. 

We note that, owing to (J2j) and r > 1, the system satisfies all the conditions in Theo¬ 
rem [273] In particular, assumption (HQ follows from (147]) . Therefore, d Xl Ui(x,t ) > 0 for 
all (x,t) G R+ x R and for all i = 1, ...,m. Now, let 

Ui,i(xi,x', t ) = Ui{x i + Z, x', t ), t ) e (— Z, cxo) x R^ 1 x R, i — 1, ..., m. 

From the boundedness of tq and parabolic estimates, letting Z —> oo upon a subsequence, 
we can assume that u it i converges uniformly on each compact set to u it00 , where (u itOC ) is 
a bounded, nonnegative classical solution of problem (JT]) in R n x R. The monotonicity 
of Ui implies that u ij0O is positive and independent of x\. We thus obtain a bounded, 
positive classical solution of problem ([T]) in R^- 1 x R. This contradicts Theorem 12.11 and 
Remark 12.11( b)) (here we assumed N > 2; the case N = 1 reduces to an ODE system for 
which nonexistence is obvious). □ 


6. Appendix 

We give the following version of the maximum principle for cooperative systems, which 
is suitable to our needs. Related results are given in [23 Section 3.8] or (TT[ Theorem 
3.2], but do not quite satisfy our requirements (unbounded domain, parabolic inequalities 
assumed on the positivity set only). Here, for given vector W := (tCi)i<j< m and real 
number M, the inequality W < M is understood as w % < M for all i — 1,..., m. 

Proposition 6.1. Let m > 2, T > 0, let Q be an arbitrary domain of R N (bounded or 
unbounded, not necessarily smooth). We denote Qt = fi x (0, T) and &pQt = (h2 x {0}) U 
(dCl x (0,T)) its parabolic boundary. Let W = (wy) G C(h2 x [0,T);R m ) D C 2,1 (Q T ; R m ) 
and denote 


Di = {(x,t) G Q t : Wi(x,t) > 0}. 
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Assume that W is a bounded, classical solution of the system 

m 

d t Wi — A uii — K\X7wi\ < Cij(x, t)wj in D t , i — 1,. .., m, (48) 

j = 1 

where K > 0 is a constant and the coefficients Cij are measurable, bounded and satisfy 

Cij > 0 for all iff- (49) 

(i) If W <0 on dpQr, then W <0 in Qt. 

(ii) Let M > 0 and assume in addition that 

m 

^^Cij < 0, i — (50) 

j = 1 

IfW<Mon dpQp, then W < M in Qp. 

Proof, (i) It follows by the Stampacchia method, e.g. along the lines of [29], Proposition 
52.21] and [29] Remark 52.11(a)]. We give the proof for the convenience of the reader and 
for completeness. 

First consider the case when 9 is bounded. Let i £ {1,... ,m}. By (148|) . we have 


[d t Wi - A Wi - K\S7wi\](wi) + < y ^Cj j (x,t)(w j )(w i ) + in Q T . 

3 = 1 

For t £ (0,T), since (wi) + (-,t) £ //](fl) by our assumption, we may integrate by parts, 
to obtain 

^4 f ( W *)+ = [ (PtWi)(.Wi) + 


2 dt 


n n TTl n 

< — / VtWi ■ V(wi) + + K / |VtUi|(iyi)+ + y^ / c ij (w j )(w i ) 
Jn Jn /=l Jn 


< 


|VK) + | 2 + /' \V(w i )+\ 2 +^- / ( Wi)l + Y^ I CpK)K)+- 

- ** j = l J 


By assumption (1491) . it follows that 

K 2 


1 d 

2 dt 


( W i) + A , / ( w i) + + / Cii(^i) + + y ^ / Cjj (iTj) + (Wi) + 


j=l 


' K 2 

< f — I - ||q- 


+ oJ2 C h(K) + + (™0 + )- 


3 = 1 


20 


Q. H. PHAN AND PH. SOUPLET 


Adding up for i — 1,..., m, we get 



for some constant L > 0. Since W < 0 at t = 0, it follows by integration that 
YJj=i fn( w i)+ < 0, hence W < 0 in Q T . 

Now, in the case of an unbounded domain, we fix £ > 0 and consider the modified 
functions 

Wi = e~ xt Wi — £0, 0 = (N + K)t + (1 + \x\ 2 ) 1 ^ 2 > 0 

with A > 0 to be chosen. We also set 

Di = {(x,t) G Qt, Wi(x,t ) > 0}. 

Since 0 t — A0 — AT|V0| > 0, we have, in D t c D,, 

d t Wi — A Wi — K\Wwi\ < e~ xt \d t Wi — A Wi — K\Vwi\ — Awy] — g[0 t — A 0 — K |V0|] 

m 

< e~ xt [-XWi + ^2 c iji x > t) w j] 

3 =1 

m m 

= [-A Wi + '^2c ij (x,t)wj\ + g[-A + y ^Cjj(x,t)]^ 

3 =1 3 = 1 

m 

< (cu(x, t ) - A)u)j + ^ c^'(x, 

J = 1 

by choosing A = max* Y1JL i 11 c *y 11 oo ■ Noting that Wi < 0 for |x| large, we may then apply 
the previous argument to get W < 0 in Q T . The conclusion follows upon letting e —* 0. 

(ii) It suffices to apply assertion (i) to the functions := Wi~ M, noting that, in view 
of assumption (l50lh we have, for each i — 1,..., m, 

d t Wi — Alii — K\Vwi\ = d t Wi — Awt — K\Vwt\ 

m m 

< ^2 dj(x, t)(wj + M ) < ^2 c ij{x, t)Wj 

3 = 1 i =1 

in D,i D {(x, t) G Qt : Wi(x,t) > 0}. □ 
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